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NOMENCLATURE 


C n (S) n times continuously differentiable on S 

p pressure 

U component velocity in x-direction 

v component velocity in y-direction 

M 1 3 2 [ ]/3x 2 + 3 2 [ ]/3 2 [ ]/3y 2 

AM ] 3“[ ]/3x 2 + 2 5!|iLl + ]/ 3y “ 

3x 3y 

3S boundary of domain S 

n. £• u 

■» implies 
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/ SUMMARY 

' 

r For more than the last 20 years there has been a concerted effort to 
solve the stationary Navier-Stokes equations; however, this has only been suc¬ 
cessful for a few special cases of primarily academic interest. An alternative 
approach has been to solve the equations numerically, and then compare the 
results with experiment. On occasion, such comparisons are in good agreement. 
However, such results are of dubious value since one has no a-priori way of 
knowing the relevance of such results until they are explicitly compared 
against experiment. Therefore, it would seem reasonable to conclude that the 
present approaches to solving the Navier-Stokes equations are of limited 
value. Accordingly, it is the purpose of this paper to show that there does, 
indeed, exist an equivalent representation of the problem that has significant 
potential in solving such problems. This is due to the fact that this equiva¬ 
lent representation of the problem consists of a sequence of Fredholm Integral 
Equations of the second kind, and the solving of this type of problem is very 
well developed. In addition, for the problem in this form, there is an excel¬ 
lent chance to also determine explicit error estimates, since one would now be 
dealing with bounded linear operators, rather than unbounded. s 


INTRODUCTION 


Ideally, one would like to obtain a numerical solution to the Stationary 
Navier-Stokes problem that is within some prescribed degree of accuracy of the 
true solution. Unfortunately, it is not possible to accomplish this objective 
by just applying some existing numerical technique directly to the Stationary 
Navier-Stokes equations, since error analysis for nonlinear equations is, for 
all practical purposes, nonexistent. However, since it is possible to replace 
the Stationary Navier-Stokes problem by an equivalent sequence of linear 
partial differential equations (ref. 1), providing some rather general condi¬ 
tions can be met, it is at least theoretically possible to obtain error esti¬ 
mates due to the vast wealth of knowledge known about linear equations. How¬ 
ever, even though attaining such error estimates is within the realm of 
possibility, one should not be misled into thinking that such a task is a 
small undertaking, for this is most certainly not the case. However, the 
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chances of someone accomplishing this task are significantly greater if the 
basic sequential problem is transformed into an equivalent form that is more 
suitable to error analysis. Accordingly, in this paper it will be proven that 
the said sequential problem can be transformed into a form that has signifi¬ 
cantly greater changes of yielding explicit error estimates. 

GENERAL DEVELOPMENT 


( 1 ) 


( 2 ) 


Mas) * t> l , M 3s > - t> 2 

* y 

Let 

vAAO + 0 AO - 0 + f, ~ f 2 

y a x y y x 

be denoted by P(0)- Hence (2) can be expressed by 

P(o) - o , 0 X (3S) « b, , ♦ OS) * b 2 

where P can be interpreted as a mapping from C I, (S) into C°(S). As demon¬ 
strated in reference 1, solving (2) is equivalent to solving the sequence of 
equations 

P(0 ) + P*(0 
n n 

3 _K 

3X dS ' 


) (il' — ) * 0 


», • vl„ • b ‘ for " ■ • • 


The Navier-Stokes equations are given by 

uu_. + vu + p - v Au + f 1 (x,y) « 0 \ 

y | in S 

uv x + w y + p y - v Av + f 2 (x,y) * 0 J 

with boundary conditions 

u(3S) « -b 2 (3S) , v(3S) = b x , 

where S is a two dimensional Green's domain with surface 3S, 

fi(x,y) e C 1 (S) , and f 2 (x,y) e C^S) , 

or equivalently, 

vAAO + 0 V ~ 0 V AO + f. - f, * 0 in S 

y a Ay x 
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providing certain conditions are satisfied (ref. 1). (i.e., primarily that 
t|» 0 , the initial guess is reasonably good), for further Insight into the sig¬ 
nificant latitude on ip 0 , the reader should consult reference 2. 


P(tf» 0 )i ] - vAA( ]+!/<„ A[ ] x + A4> 6 f ] - A* 0 ( ] x - A( ] 


Equivalent to equation (3) is the problem 


P <V + p, <*n>*n-° * TT 
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• 0 , n * 0,1,2, 
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by letting 
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Now consider (4) under a slight change in boundary conditions, in particular: 


*(♦_) + P* (*„)** ' 0 

m mm 
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3n 


as 


oi * 0|1|2) • • • 


(5) 


If is a solution to (5), it directly follows that is also a solution 

of (4), and vice versa (i.e., equation (4) ■* y m | * constant, without loss of 

I 3S 

generality assume constant * 0). Therefore, it is sufficient to focus our 
attention on equation (5), which in detail is given by 


uAAiJ/ + 
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- * m Mm + P(* m > • 0 
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With 


* n OS) * 0 




0 , m - 0,1,2, . . . 


Hence, 

"+m * v J[ G( \ A \ + ' \ A \ )ds +if s PG ds (6) 

where G is the Green’s function of the biharmonic Equation [3] 

AA$ * 0 in S, with $(3S) * 0 , $ (3S) ■ 0 

n 

Lemma: 


Proof.O 


(6) - (5) 









For convenience and clarity, equation (6) will be expressed in the form: 

$ • A(4i n )^ + f<* m ) (6) 

where 

*<*„>' 1 - i f 0(* A| ] x t 4*^1 ) y - A*. I )„ - *. A[ ) y )d» 

and 

f(* m ) - - J PU>JG ds 
m v 's m 

Lemma: If the hypothesis of theorem 1 of reference 1 is satisfied, then there 
exists a solution of equation (6) for all m and 


n 

4>* - Urn £ - * # 

n-»® i*o 

where ij<* is the unique of equation (2). 

Theorem: Under the condition of the above lemma, is a solution to the 
Fredholm Integral Equation of the second kind 


^*(x',y') + J ^(x.yJK^x* ,y‘ ,x,y)dx dy «= fU^Hx'.y') 
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where 
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-AGy(x',y' 


» x >yHm < x *y> + 4G xV 

x y 


+ 



+ ^m xy (G xx - G yy ) 

Proof: Since will be a solution of (6) under the hypothesis, it is suffi¬ 

cient to show that a solution to (6) is alos a solution of (7). For clarity, 
denote (6) in the form 


- A(^ 0 H* *= f(<{> 0 ) 


or equivalently. 
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By Green’s Lemma 
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Therefore 
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4 .* + — 

V 

/J[4- y )J ; * + c KV* ■ *KV‘ 

- (° *%) 

4 1 
X 1 

4 ,* 

l+ i/ s »*K c \) x + 

(g Aip g ) - A(gv 0 ) 

' X'y ' X 'y 

- ( G 4 »*y) 

JdS 

x-* 

K W "y) x 

+ ( G W’o ) - A ( G * 0 ) 

' X'y ' X'y 

- ( G a \)J 



- “(v. 

+ G « 0 \ + G A 4 0 + 

j xy/ 7 x 

G Ai}, g - a(g ip 0 + 

xy ' 1 x 

W. ) - C X 

xy' 

A <' 0 

• 4x<\ 

, ■ V° x ) + ( Gy A< ’°* " 

G X *%) 




dS •= f( V 0 ) 


- f(v 0 ) 


- G A <[ 0 
y xy 


4 (vo) - (mo) + ( G x n ) 

v y' v y'xx ' y'yy 

■ (v V + =x* v ) x - ( G xy*. y + V«J y 

• G XXX*t y + G * 4 » y> , y + C X*»xxy + G *»y*'y + 2 ( GxX '"xy + G »y l| '»yy) 


5 




- a (vO ■ -(v.,) xx - (s*. x )„ 
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Therefore, (6) reduces to 
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Hence, (6) ** (7). 

CONCLUSION 

Therefore, under the conditions as cited, the sequence of solutions 
generated by (7) converge to .the solution of (1) in the following sense 

"t », 

i°o 

fl * if>* - n, 

o A 

Therefore, 
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Hence the original problem, equation (1), which falls within the frame¬ 
work of nonlinear operator theory, an area that little is known, has been 
replaced by equation (5), which falls within the framework of bounded linear 
operator theory, an area for which there exists a vast wealth of information. 

Of course, it could be argued that it is theoretically possible that the 
sequential representation of the solution converges so slowly that the results 
are of questionable value. However, from all indications this will not be a 
problem, for as pointed out in reference 1, there are parameters at one's 
disposal in the method that can be adjusted to speed up the rate of conver¬ 
gence; in fact it was demonstrated in reference 2 that by judicious selection 
of the variable parameters available they were able to get numerically adequate 
convergence with just a few iterations. 
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